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Majorana fermions were proposed to occur at edges and interfaces of gapped one-dimensional 
systems where phases with different topological character meet due to an interplay of spin-orbit 
coupling, proximity-induced superconductivity and external magnetic fields. Here we investigate the 
effect of strong particle interactions, and show that the helical liquid offers a mechanism that protects 
the very existence of Majorana edge states: whereas moderate interactions close the proximity gap 
which supports the edge states, in helical liquids the gap re-opens due to two-particle processes. 
However, gapless fermionic excitations occur at spatial proximity to the Majorana states at interfaces 
and may jeopardize their long term Majorana coherence. 
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I. INTRODUCTION 

Majorana fermions (MFs) were introduced in 1937 as 
neutral particles that are their own antiparticles, with the 
aim of describing neutrinos.^ More recently an increasing 
number of candidate MFs were suggested as quasiparti- 
cles in condensed matter systems, including the Moore- 
Read quantum Hall state^ at filling fraction 5/2, which 
can be described as a p wave superconductor of compos- 
ite fermions.^ In this picture MFs are zero energy bound 
states localized at vortices^ and show non-abelian statis- 
tics upon braiding those vortices;^ for a recent review see 
Refs.|6|and|7l 

Searching for alternative experimentally feasible can- 
didates, inspired by the one-dimensional (ID) model of 
Kitaev,!^ several groups have proposed physical realiza- 
tions of MFs as edge states of ID systems, including 
electrostatic defect lines in superconductors,'^ semicon- 
ductor quantum wires proximity coupled to a supercon- 
ductor ,^^21111 q^agi. ID superconductors,^^ and cold atoms 
trapped in ID.^^ In semiconductor wires these edge states 
can be controlled by tuning external gates, and networks 
of such wires are e nvisio ned to perform non-abelian quan- 
tum computationP^^ 

The practicability of such ID applications critically 
hinges on the stability of Majorana fermion states against 
particle interactions. Focusing on the most elementary 
realization of a quantum wire with a single spin polar- 
ized fermion band proximity coupled to a superconduc- 
tor and with spin-orbit coupling, which was showii^ to 
reduce to Kitaev's modeP - Gangadharaiah et al.^ in- 
deed have shown that even moderately strong interaction 
may compromise the stability of Majorana fermions: be- 
yond a certain strength, interactions remove the proxim- 
ity gap, and along with it any accompanying Majorana 
bound states. Here we consider the impact of interactions 
on the full 'helical liquid', i.e. a system of two counter- 
propagating fermion bands carrying opposite spin, as re- 
alized as the surface state of two-dimensional topological 
insulators,^'' or proximity coupled semiconductor quan- 
tum wires subject to spin-orbit interaction,!^^ cf. Fig.[l] 
We show that Majorana fermion states forming in these 



systems enjoy a much higher degree of stability. Concep- 
tually, the resilience of the helical liquid to interactions is 
rooted in momentum conserving two-particle scattering 
processes between the constituting fermion bands .^J*^ 
Such processes open gaps that resist interactions and ul- 
timately lead to the stabilization of Majorana states. At 
the same time, strong interactions may act as a source of 
Majorana fermion decoherence. 
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Figure 1. Helical liquid formed by (a) surface bands of a 2D 
topological insulator, or (b) the low-momentum excitations of 
a quantum wire subject Rashba interaction. Inset: Gaps are 
induced both by magnetic field (B) or superconductivity (A). 
While for /j, < \/B'^ — A^, A < B and no interactions, the gap 
is controlled by the superconducting gap dominates for 
/i > ^JB^ — A^ and a different topological sector is obtained. 



The outline of the paper is as follows. Having intro- 
duced the generic model for an interacting helical liquid 
in Sec.|TI|and mapped it into a spin chain model, we map 
out the helical liquid's global phase diagram in Sec. |III| 
in a parameter space spanned by proximity coupling, ex- 
ternal magnetic field strength, interaction strength and 
chemical potential. At strong interactions the system 
turns out to support a phase of gapless fermionic exci- 
tations in close parametric neighborhood of the Majo- 
rana/Ising quantum critical point. In a quantum wire 
subject to slowly changing parameter profiles, considered 
in Sec. IIVI this means the existence of a domain of low 
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lying excitations in close spatial proximity to a localized 
Majorana state. As pointed out in the concluding section 
\V\ such excitations may spoil the long term coherence of 
the Majorana fermion, however fermion parity protected 
readout schemes may b e emp loyed to probe a certain de- 
coherence free subspace.!^^^ Details on the derivation of 
the phase diagram using bosonization methods are in- 
cluded in the appendix. 



II. MODEL 

A helical liquid is defined by the Hamiltonian Hq = 

^0 i^Lii^Fidx - lj)i^Li + il^^^^i-vpidx - lj)i^R^, (1) 

where ^pR^{^l^Li) are right (left) moving spin up (down) 
fermion fields, vp is the Fermi velocity, and /i the chemi- 
cal potential. Gaps of different topological signature may 
be realized by coupling the system to a magnetic field B 
(Zeeman gap), and to an s-wave superconductor (prox- 
imity gap), 

5n = Bi^l^i^R^ + A^Li^R^ + H.c; (2) 

see Fig. 1. To illustrate the nature of the transition 
between the gapped phases, we consider for simplicity 
the limiting case ja = and introduce two Majorana 
fields tpLiix) = {ixi{x)^X2{x))/V2, ^l^n^ix) = {xi{x) + 
iX2{x)) / V^, in terms of which the Hamiltonian reads 

no^6n = Y.]^i {Xji^d^Xj - Xji^d^Xj + ^rUjXjXj). 
with mi ^2 = A ^ B (with gauge choice A > 0). At 
\B\ = A, one of the two Majorana modes becomes mass- 
less which is a signature of a bulk Majorana/Ising quan- 
tum phase transition (QPT).^^ Tuning parameters in real 
space, so that mi{x) changes sign, one obtains a Majo- 
rana zero mode localized around xq with mi(xo) = 0, 
as can be verified by an explicit solution of the inho- 
mogeneous quadratic theory. The Majorana edge modes 
forming at the interface between distinct phases are topo- 
logically protected and, hence, candidate quantum bits. 

Here we will focus on the impact of interactions on 
the QPT, B A. To see that sufficiently weak interac- 
tions have no qualitative effect, we consider the Majorana 
representation of an interaction term, Hint ^ X1X1X2X2, 
and treat the gapped Majorana component as a c-number 
X2X2 {X2X2)' The ensuing mean field approximation 
can be absorbed into a redefined mass parameter mi of 
the critical mode, i.e. a weakly shifted transition point. 



A. Generic interactions and mapping to a 
spin-chain model 

The helical liquid permits two types of time reversal in- 
variant interactions, forward and umklapp (two-particle) 



scatterin^I^ 

^um ^nV^l;9xV^l;V^i?t^^V^i?t + H.C. (3) 

To make progress with the interacting model, we map the 
Hamiltonian = Ho + Hfw + '^um + to the XYZ spin 
chain model with both non-staggered and staggered mag- 
netic fields, HxYZ = ^i^i (up to a constant), where 

ni= Yl JaStsf^,-[^i + B{-iy]s!, (4) 

a=x,y,z 

and the coupling constants J^^y = J±A>0, J = vp, 
and Jz > are fixed by the condition that after a Jordan- 
Wigner re-fermionizatiorl^ 

and expansion of the lattice fermions aj in terms of the 
left and right movers, 

aj - e'^'^iljR^ix) + e"^^^V^Li(^), 

(with X = aj and lattice constant a ^ 1), one recovers 
the starting Hamiltonian H with Jz = ^ = ^ = Su- 
We emphasize that the details introduced by the lattice, 
such as the above relation between the interaction cou- 
plings g2 , g4 and gu , are unimportant for an identification 
of the morphology of the global phase diagram and the 
universality of the QPTs;^^ on the other hand commen- 
surability with the half filled lattice at /i = should not 
be disregarded as a lattice artifact, since it refiects the 
physical umklapp process which involves the original zero 
momentum helical modes; see Fig. l."^ 



III. PHASE DIAGRAM 

To start our discussion of topologically distinct phases, 
we consider the phase diagram spanned by the parame- 
ters, J^, 5, and A, starting with the case /i = 0, cf. 
Fig. [2] In the non-interacting limit, = 0, and at the 
particular parameter configuration A/ J = 1, the Hamil- 
tonian Q is seen to reduce to a transverse field Ising 
model (after a tt rotation of each second spin around x), 

^ 2A5f - BSt. (5) 

The latter supports an Ising type transition at 1^51 = A 
(blue diagonal lines in the = plane of Fig. [2|: at 
\B\ < A the Ising symmetry —S^ is spontaneously 

broken yielding a doubly degenerate Neel-x(2)-ordered 
phase, where (n) denotes a phase with degeneracy n. At 
\B\ > A di unique Neel-z(l) state is chosen by the direc- 
tion of the staggered field B. The phase transition is long 
known to have a Majorana critical mode,^^ as identified 
above along the \B\ = A line. In line with the mean 
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field argument above, finite > is irrelevant at this 
Majorana QPT. 

To understand what happens at generic interaction 
strength, we first consider the line A = B = (horizontal 
red line in Fig. [2| , where the model reduces to an XXZ 
spin chain with spin-z rotation invariance. This model 
contains an XXX Heisenberg critical point separating a 
phase of gapless planar spin fluctuations at < J [cor- 
responding to a Luttinger liquid (LL) in the fermion rep- 
resentation] from a doubly degenerate Neel-z(2)-ordered 
phase at > J; expressed in terms of the universal 
Luttinger interaction parameter, the transition occurs at 
K — 1/2, where K = 1— to first order in the param- 
eters ^2,4,w'^In fermionic language the order parameter 
in the Neel-z(2) phase reads O = ipl^^^ipR^ -\- H.c, i.e. fer- 
romagnetic ordei^^ (of the physical spin) in a direction 
parallel to the external field. 




interactions [Jz] 
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relevant. Specifically, finite A perturbes the Hamilto- 
nian by SHi = SfSf_^i — Sf S^_^i, which breaks spin ro- 
tation invariance and carries positive scaling dimension 
xa = 2 — 1/K. This operator causes flow towards the 
doubly degenerate Neel-x(2) order. Finite B couples to 
the spin chain through a staggered magnetic field in z- 
direction, leading to an operator with positive scaling di- 
mension xb = 2 — and drives the system towards the 
non-degenerate Neel-z(l) phase. The discussion above 
implies that the line of QPTs A = 5 at = evolves 
into a 'surface' of transition points (blue surfaces ema- 
nating from the segment Jz/ J < 1 with l/2<i^<lon 
the Jz axis in Fig. [2|. The cusps of these surfaces for 
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follow from the difference in scaling dimensions of the 
two competing perturbations, xb > xa ^oi K < 1. 

At strong interactions, K < 1/2, weak A and B be- 
come irrelevant perturbations and cease to affect the 
Neel-z(2) phase. To understand what happens in this re- 
gion, it is important to notice that the Heisenberg point 
i^ = l/2,5 = A = 0is terminal to a second line of XXZ- 
models, specified by J^^ = or A = — J > 0, 5 = 
(diagonal red line in the top panel of Fig. [2|. This line 
corresponds to gapless planar fluctuations around a con- 
served spin-y, and its critical properties are equivalent to 
those of the horizontal line save for an exchange z ^ y 
of the invariant axis.l^ As follows from a bosonization 
analysis of this diagonal line carried out in the appendix, 
the perturbation A is again relevant with respect to the 
LL line A = — J > 0, driving the system either to 
the Neel-x(2) phase for A > J;^ — J, or to the double 
degenerate Neel-2:(2) ordered phase, for A < — J; the 
staggered field B being now perpendicular to the con- 
served spin-axis is also relevant and drives the system to 
the Neel-2:(1) phase. We infer that the Majorana transi- 
tion sheets extrapolate to strong interactions and merge 
with this LL as shown in the figure. 



Figure 2. Top: schematic phase diagram of the XYZ model, 
Eq. Q, at /i = 0. Neel ordered phases are denoted by a(n), 
where a = x,y,z denotes the magnetization axis and n the 
ground state degeneracy. At the thick (red) lines the model 
reduces to a critical XXZ chain described by a LL theory. The 
blue surfaces are Ising transitions described by MF bulk crit- 
ical modes, between the doubly degenerate Neel x(2) phase 
and the nondegenerate Neel z(l) phases. The latter are con- 
nected with the doubly degenerate Neel z{2) state occurring 
at B = 0, J2 > J+ |A|. Bottom: schematic phase diagram at 
fi ^ 0. Inset: phase diagram of the XXZ model at A — B — 0. 
In the dark shaded region, where 1/2 > K > 1/4, the LL 
phase is stable against small A. 

Perturbations in both finite B and A around the weak 
interaction segment, 1 > K > 1/2, are analyzed in 
the appendix using bosonization and are found to be 



A. Phase diagram, (/i 7^ 0) 

At finite /i the horizontal line 5 = A = of gapless 
excitations at < J opens to become a surface in the 
A = plane. This is easily understood in fermionic lan- 
guage where it means that for /i ^ a finite magnetic 
field is needed to create an excitation gap in the heli- 
cal liquid as is apparent from the inset of Fig. 1. For 
larger interactions the threshold field becomes smaller, 
which can be seen as a precursor phenomenon of ferro- 
magnetism. The behavior at the threshold has a uni- 
versal description in terms of commensurate to incom- 
mensurate (C-IC) transitions; see the appendix for more 
details. At values of B large enough to sustain a gap, 
a surface of transition points Acnt = f{B^Jz) separates 
phases of Neel z and x order, as qualitatively shown in 



4 



Fig. [2] (blueish surface.) Importantly, the system sup- 
ports a second domain of gapless excitations, shown as a 
tubular structure in Fig. |2j 

To understand this region, one has to notice that as 
one sweeps the interactions in the entire critical region 
at 5 = A = (bolded line with an end arrow), where 
Eq. Q becomes an exactly solvable XXZ chain in a lon- 
gitudinal magnetic field, it is known that K drops from 
1 to 1/4 at a C-IC transition (as opposed to K = 1/2 at 
the Kosterlitz-Thouless transition obtained for /i = 0, cf. 
inset of Fig. |2j) This fact entails the existence of a finite 
interval of values for which the LL is strongly repul- 
sive 1/2 > K > 1/4 such that A is irrelevant, xa < 0, 
i.e., there exists a three-dimensional region of parameter 
values as indicated in the figure by the red tube for which 
the system remains gapless, and equivalent to a Luttinger 
liquid in the language of fermions. As we argue in the ap- 
pendix, this critical region is expected to extend to large 
values of A. Whereas such an interaction induced gap- 
less phase was predicted in the spin polarized quantum 
wires, in the helical liquid considered here upon fur- 
ther increasing interactions the umklapp term becomes 
relevant at the C-IC transition at the curve shown as the 
diagonal red line in the bottom panel in Fig. [2j and the 
gap re-opens. In the appendix we demonstrate that this 
C-IC transition line connects with the Majorana transi- 
tion sheets. 



IV. LOCALIZED STATES AT INTERFACES 
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Figure 3. Spatial variation of system parameters at weak (a), 
or strong (b) interactions. Variation profiles correspond to 
paths a ^ h ^ a 01 a' ^ h' ^ a' oi Fig. [2] respectively. 
The doubly degenerate x(2) center region is associated to 2 
unpaired MFs (filled red squares) located near the interfaces. 
For strong interactions spin fluctuating regions form near the 
interfaces, which in fermion language corresponds to LLs. The 
ensuing low energy modes are coupled to the MFs. 



Having mapped out the phase diagram, we next discuss 
the implications for the Majorana edge states. To this 
end, imagine a space dependent change of parameters, 
e.g., along the path a ^ h ^ a m Fig. [2] Within the 
spin picture, this leads to a Neel z(l) x{2) z{l) 
ground state structure, cf. Fig. |3] upper panel. The 
central x{2) domain breaking the Ising symmetry 
—S^ is two-fold degenerate, and at the same time lacks 
any local order in fermion language. This signifies the 
presence of unpaired zero energy MFs at the interdaces. 
This correspondence can be seen by inspection of the 
Jordan- Wigner lattice fermion system corresponding to 
the center Neel x{2) region (sites 1 — N). Introducing 
a pair of Majorana fermions for each lattice fermion as 

lA,j = CLj + CLj^ and jbj = — ^(%' ~ ^])^ finds that 
the two terminal Majorana fermions do not enter the 
Hamiltonian.^ The coupling between those MFs, 

with Pauli matrices ctj = 25^, is an operator that fiips 
an entire Neel x domain, and therefore is suppressed ex- 
ponentially in N. 

More interesting things happen along a sweep of pa- 
rameters in the strongly interacting system. Consider, 
for example, the path a' ^ h' ^ a' in Fig. [2j where the 



proximity coupling, A is increased in a region of space 
to pass from a region with an interaction gap {a') to 
one with a proximity gap {h'). This corresponds to an 
z(2) x(2) z{2) ground state structure (Fig. [s] bot- 
tom) where, however, the z{2) degeneracy of the outside 
regions does not imply extra Majorana states because it 
corresponds to the degeneracy due to local ferromagnetic 
order in terms of electrons. This shows that the very 
existence of MFs at interfaces is robust against strong 
interactions. However, upon moving from a' through the 
phase boundary towards h' ^ one crosses the tubular crit- 
ical region of Fig. [2] In spin language, this is a region of 
gapless planar fiuct nations, in fermion language a LL. Ei- 
ther interpretation shows that in close spatial proximity 
to the Majorana fermion state a (Luttinger) liquid of low 
energy excitations forms. The spatial extension of this 
region, Lll, is the larger the more shallow the parame- 
ter profile of the wire is. While the velocity of excitations 
vanishes at the C-IC transition, we nevertheless find by 
solving an effective Schrodinger equation, that the finite 
size gap scales as ^ 'i^eff/^LL, where '^eff = Oivp)- In 
comparison, the energy gap towards higher fermionic ex- 
citations f or weak interactions is much larger and scales 
aJ^ ^ ^J ttiVq^ I L where m/L is the slope of the gap 
and Lll ^ L. Whether or not fiuctuations of the LL 
get excited primarily depends on the ratio T/(i'eff/^LL) 
between excitation energies and temperature, T. 
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DISCUSSION AND SUMMARY 



s:iven m Eqs. (1,2,3) reads^a 



What is the impact of the LL on the nearby bound- 
ary Majorana states? Unhke non-fermionic low energy 
quantum fluctuations (phonons, nuclear spins, etc.), the 
presence of the LL may spoil the long term quantum co- 
herence of Majorana fermion states. To see this, denote 
the MF operator of a phase crossing by 7^. Non-fermionic 
excitations do not change the parity of the overall MF 
number, and in the absence of low energy fermion states, 
they can couple to MFs only via Majorana bilinears ijijj. 
In view of the large separation between Majorana sites, 
such operators are exponentially suppressed. By con- 
trast, the coupling to a system of fermionic low energy 
excitations is mediated by operators of the form 7^(0+0'*'), 
where is a fermion creation operator. This causes en- 
tanglement of individual MFs with a collective quantum 
environment, and may compromise long term quantum 
coherence. It may become an issue at strong interac- 
tions in quantum operations based on the tunnel cou- 
pling of endpoint MFs in systems of connected quantum 
wires .SS' However, the conservation of the parity of the 
total number of fermions in the interface allows to define 
a protected Majorana subspace,^^ which can be probed 
by more elaborate interference measurement.'^ 

Summarizing, we have mapped out the phase diagram 
of interacting helical liquids subject to both proximity 
coupling to a superconductor and a magnetic field. We 
have reached the conclusion that even in parameter re- 
gions where the interactions are strong enough to open 
a gap on their own account, the increase of a compet- 
ing proximity gap will generate Majorana excitations. In 
general, however, these excitations will suffer decoher- 
ence, and this is due to the fact that in close proximity 
to the phase boundaries between interaction and prox- 
imity gaps, the system supports low energy electronic 
excitations. It will be interesting to confirm this picture 
numerically. 
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Appendix A: Derivation of the phase diagram from 
bosonization 



In this appendix we use bosonization to reproduce the 
phase diagram (PD) in Fig. 2. The bosonized version of 
the helical liquid Hamiltonian H = 'Ho +'^fw +'Wum + ^^'^ 



(Al) 



— -^^ cos(Vl67r^) + — cos{^/47^(j)) — — cos(A/47r6>), 

where v = '^f + + ^(^i^^l^^n) renormalized 
velocity, K is the LL parameter given above to first or- 
der in g2^g4.-,gui the Bose fields satisfy [(j){x)^dx'0{x')] = 
i5{x — x')^ and the helical electrons are represented by 
i^RULi) ^ is ^iso contin- 

uum theory of the XYZ spin chain Eq. (4), where spins 
are represented b}Eil 



S+ix) 



-^dx(l){x) + ^ cos[>/47r0(x) — TTx], 



27r 



{(-!)" + cos[V4^(/)(x)]}. (A2) 



1. Commensurate case (/x = 0) 

At 5 = A = the standard renormalization group 
(RG) analysiJ^ shows that for infinitesimal gu = 0+ the 
umklapp term becomes relevant upon increasing interac- 
tions beyond a critical value. This occurs at a Kosterlitz- 
Thouless (KT) transition when K = 1/2, and the system 
goes from a gapless LL at > 1 /2 to a gapped phase at 

< 1/2 with a pinned value of the bosonic field (0) = 
or ((/)) = y^/2, which classically minimize the gu term 
in Eq. (TaTI). From Eq. (|A2|), this phase has Neel order 



along and corresponds to the Neel z{2) phase in the 
PD at large interaction. 

The perturbations A and B around the weak inter- 
action LL segment, 1 > K > 1/2, are both relevant.!^ 
Specifically, infinitesimal A carries positive scaling di- 
mension X/\ = 2 — 1/K and leads to pinning of the 
field at (0) = or (6) = y^, minimizing the correspond- 
ing term in Eq. (|A1|). Using Eq. (A2), this phase has 



Neel-order along x. Analogously, infinitesimal B > 
leads to an operator with positive scaling dimension 
xb = 2 — and drives the system towards a phase 
with ((/)) = \pKl2 (or ((/)) = for 5 < 0), associated 
with Neel-order along z. Those two Neel z states driven 
by 5 > or 5 < are seen to correspond to the two 
Neel 2: (2) states driven by strong interactions. This im- 
plies that the Neel z(X) and 2: (2) phases are connected 
to each other in the PD without crossing any QPT. This 
RG picture valid around the LL segment 1 > > 1/2 
(i.e. < J < in terms of spin chain parameters) im- 
plies that the Ising surfaces inferred to emanate from the 
noninteracting = plane at \B\ = A, merge with the 
LL segment 1 > > 1/2 at 5, A ^ 0. 

At strong interactions there is another critical LL line 
at A = — J > depicted as the diagonal red line in the 
top panel in the PD, described by an XXZ model with 
Jx — Jy (similarly, there is a third LL line at A = 
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J — < with Jy = Jz > Jx, not shown). The LL the- 
cry is given by 'Hj^=j^>j^ = \ {^{d,<P'? + K'{dJ'f) 
where t he b osonic fields (j)' and 0' are related to spins 
by Eq. ( A2 ) but with conserved axis being now y rather 



than namely 

sy < 



/tt 



27r 



27r 
[(-1)^ 



cos(v47r0' — TTX), 

■ + cos(V^(^0]- (A3) 



Here, the LL parameter K' and velocity v' are given \y^^ 
K' = 



2(7r- 



and V 



2 arccos (Jy/ J^c) 



Different 



than the J < JzLL segment, here the staggered field B is 
perpendicular to the conserved axis, and using Eq. (A3) 



we find that the most relevant perturbations from this 
LL line are 



/o A / 

6H = cos{^/ne') + — cos(V^i90, 

^/7^ TT 



(A4) 



where A' = A— 7^+ J. The second term, as the term cx A 
in Eq. (Al), is the leading U(l) breaking perturbation 
caused by finite A' which renders Jx ^ J^. T he scaling 
dimensions of the two operators in Eq. (A4) are x^^ = 
2 - l/(4i^0 and x'. = 2-l/K', respectively. Both 
operators in Eq. ( A4 ) are relevant long the /S. = Jz — J > 



line, tending to localize the 0' field. They may or may 
not compete with each other depending on their relative 



or 



sign. Positive A' locks 6' to {6') 
corresponding to the Neel x{2) phase, whereas negative 
A' locks 0' to {0') = or {0') = y^, corresponding to 
the Neel z{2) phase. Upon increasing B at fixed A' < 
[starting from the z{2) phase below the diagonal red 
line PD] there is no transition since the minimas of both 
cosines in Eq. ( |A4[ ) overlap. On the other hand, starting 
at fixed A' > [at the x{2) phase above the diagonal 
red line in the PD] there is a transition upon increasing 
B to either {0') = for 5 < or {6') = for 5 > 0, 
corresponding to t he N eel z{l) states. Such transition 
has Ising character .'^^'^ This RG picture implies that the 
Ising surfaces discussed above should eventually merge 
with the LL line A = Jz — J > 0. Close to this line 
the competing operators with scaling dimensions x'^ > 
x'^ lead to the cusp singularity A^^^^ oc 5^a/^b which 
connects to Eq. (IgI) at the Heisenberg point = Jy = Jz- 



2. Incommensurate case (/i ^ 0) 

We initially consider the A = plane, and fix B and 
Jz. Large enough /i, acting as a uniform magnetic field 
on the spin chain, w ill le ad to a finite magnetization, 
(S^) ^ 0. From Eq. (A2), this means that the phase (j) 
will acquire a linear term (j){x) = y/7r{S^)x + where 
... denote fluctuations, therefore the B and gu terms in 
Eq. (Al) become oscillating in space and can not open 



a gap. Thus the system is in an incommensurate LL 
phase in which {S^) changes continuously as function of 
/i. Upon decreasing jn the magnetization decreases till 
(S^) ^ 0, at which point the B and gu terms no longer 
oscillate, and the system enters into either one of the 
commensurate z{l) or z{2) phase for either 5 ^ 0, or for 
B = and sufficiently large Jz , respectively. This follows 
from the stability of the gapped phases in the /i = 
PD against finite fi. Hence there is a commensurate to 
incommensurate (C-IC) transition between the z{l) or 
z{2) phases and the LL gapless state as function of /i. 
Alternatively, at fixed /i this transition occurs as function 
of B or Jz as shown in the PD. The value of the LL 
parameter at a C-IC transition is renormalized to the 
universal value K = where n is the degeneracy of 

the commensurate phase. Hence, we infer the existence 
a critical LL surface in the A = plane in the PD at 
/i ^ 0, with LL parameter K = 1 dit its edges, except at 
the B = tip, where K = 1/4. 




Figu re 4 . Schematic phase diagram of the ANNNI model 
Eq. (A5). Sarting at the Neel-x(2) phase, upon increasing the 
interaction Jz at fixed A//i (trajectory along the dashed line) 
one first crosses a KT transition into a gapless LL phase, and 
then one crosses a C-IC transition into the Neel-z(2) phase. 
This trajectory corresponds to crossing the red tube in Fig. 
2 at A = J, B = 0, /X ^ 0, upon increasing J^. The ANNNI 
model also has an Ising (I) transiton between the Neel-x(2) 
phase and a polarized phase. 



Using the field theory Eq. (Al) describing the critical 
surface at the A = plane, we may analyze the effect of 
an infinitesimal A. Since xa = 2 — l/iC, this perturba- 
tion is relevant for K > 1/2, and in this case the field 
is pinned and the system enters the Neel-x(2) phase for 
infinitesimal A. However A is irrelevant for K < 1/2. 
Since K varies continuously in the LL surface at A = 
except at the the 5 = tip with K = 1/4, we con- 
clude that there must exist a finite critical subarea with 
K < 1/2, as shown schematically by the red dark area 
in the A = plane in the bottom panel of the PD. This 
critical subarea is stable against A and hence extends 
into a three dimensional critical region above the A = 
plane up to some finite value of A, as shown schemati- 
cally by the tube in the PD; as follows from the universal 
values of the LL parameter at C-IC transitions quoted 
above, we conclude that K = 1/2 at the 2D boundary 
of this tube with the Neel x{2) phase, and K = 1/4 at 
its C-IC transition with the Neel z{2) phase. At 5 = 0, 
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the C-IC transition is the only region in the PD which 
is sensitive to smah B. Hence, the Majorana transition 
surfaces must merge with this C-IC transition hue. 

It is natural to expect that this 3D manifold extends to 
large values of A along the critical line A = — J (for 
/i ^ the critical manifold will shrink to this critical 
line). As an evidence for this expectation, we observe 
that for 5 = and A = J Eq. 4 reduces to the extensively 
studiec^ ANNNI (Axial Next Nearest Neighbor Ising) 
model, 

^ 2ASfS'f+, + J.StSt^^ - iJiSt, (A5) 



whose phase diagram is shown schematically in Fig. [4] 
The so-called "floating phase'" of this model (red shaded 
area), can be identified with the LL phase inside the criti- 
cal tube in Fig. 2 at finite /i. Note, however, that whereas 
there exists an analytic proof that this LL phase emerges 
out of the A = 0, /i = multicritical point in Fig.[4p^^ 
there is still no full evidence that this LL phas e extends 
to the x{2)-z{2) transition point at /i = O.EIEll 



^ E. Majorana, Nuovo Cimento 14, 171 (1937). 
2 G. Moore and N. Read, Nucl. Phys. B 360, 362 (1991). 
^ N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 
^ N.B. Kopnin and M.M. Salomaa, Phys. Rev. B 44, 9667 
(1991). 

^ G.E. Volovik, JETP Lett. 70, 609 (1999). 

^ F. Wilczek, Majorana returns, Nature Physics 5, 614 

(2009) . 

^ C. Nayak, S. H. Simon, A. Stern, M. Freedman, S. Das 

Sarma, Rev. Mod. Phys. 80, 1083 (2008). 
^ A. Y. Kitaev, Physics-Uspekhi 44, 131 (2001). 
^ M. Wimmer, A. R. Akhmerov, M.V. Medvedyeva, J. 

Tworzydlo, and C.W. J. Beenakker, Phys. Rev. Lett. 105, 

046803 (2010). 

R. M. Lutchyn, J. D. Sau, and S. Das Sarma,Phys. Rev. 
Lett. 105, 077001 (2010). 

Y. Oreg, G. Refael and F. von Oppen, Phys. Rev. Lett. 
105, 177002 (2010). 
^2 A. C. Potter and P. A. Lee, Phys. Rev. Lett. 105, 227003 

(2010) . 

L. Jiang et al, Phys. Rev. Lett. 106, 220402 (2011). 

J. Alicea, Y. Oreg, G. Refael, F. von Oppen and M. P. A. 

Fisher, Nature Physics, doi:10.1038/nphysl915. 

J. D. Sau, S. Tewari and S. D. Sarma, Phys. Rev. A 82, 

052322,(2010). 

S. Gangadharaiah, B. Braunecker, P. Simon and D. Loss, 

Phys. Rev. Lett. 107, 036801 (2011). 

M. Konig et al, Science 318, 766 (2007). 

In contrast, in the polarized regime studied in Ref.lTBltwo- 



particle processes require special unrealistic filling. 

C. Xu and J. E. Moore, Phys. Rev. B 73, 045322 (2006). 

C. Wu, B.A. Bernevig, and S.-C. Zhang, Phys. Rev. Lett. 
96, 106401 (2006). 

J. B. Zuber and C. Itzykson, Phys. Rev. D 15, 2875 (1975). 

22 A. R. Akhmerov, Phys. Rev. B 82, 020509 (2010). 

23 F. Hassler, A. R. Akhmerov, C.-Y. Hou and C. W. J. 
Beenakker, New J. Phys. 12, 125002 (2010); E. Grosfeld, 
B. Seradjeh and S. Vishveshwara, arXiv:1004.2295 (unpub- 
lished) . 

2^ T. Giamarchi, Quantum Physics in One Dimension (Ox- 
ford University Press, New York, 2004). 

2^ In quantum wires Qu is typically smaller than ^2,4 as it is 
generated by a combination of interaction and spin-orbit 
coupling. 

2^ H. J. Schulz, Phys. Rev. B 34, 6372 (1986). 
2^ G. Delfino and G. Mussardo, Nucl. Phys. B 516, 675 
(1998). 

2^ For a review of the ANNNI model see, e.g., W. Selke, Phys. 

Rep. 170, 213 (1988); B. K. Chakrabarti, A. Dutta and P. 

Sen Quantum Ising Phases and Transitions in Transverse 

Ising Models (Berlin: Springer, 1996). 
29 J. Villain and P. Bak, J. Phys. (Paris) 42, 657 (1981). 

P. Rujan, Phys. Rev. B 24, 6620 (1981); see also E. Sela 

and R. G. Pereira, Phys. Rev. B 84, 014407 (2011). 

D. Allen, P. Azaria and P. Lecheminant, J. Phys. A: Math. 
Gen. 34, L305 (2001). 

^2 M. Beccaria, M. Campostrini and A. Feo, Phys. Rev. B 
76, 094410 (2007). 



